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Abstract 

We investigate the validity of the generalized second law (GSL) of gravitational thermo- 
dynamics in the framework of f(T) modified teleparallel gravity. We consider a spatially 
flat FRW universe containing only the pressureless matter. The boundary of the universe 
is assumed to be enclosed by the Hubble horizon. For two viable f(T) models containing 
f(T) = T + fi 1 (-T) n and f(T) = T-/x 2 T(l-e^^), we first calculate the effective equation 
of state and deceleration parameters. Then, we investigate the null and strong energy con- 
ditions and conclude that a sudden future singularity appears in both models. Furthermore, 
using a cosmographic analysis we check the viability of two models. Finally, we examine the 
validity of the GSL and find that for both models it is satisfied from the early times to the 
present epoch. But in the future, the GSL is violated for the special ranges of the torsion 
scalar T. 
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1 Introduction 



Observational data coming from the type la supernovae (SNela) surveys, large scale structure 
(LSS), and cosmic microwave background (CMB) anisotropy spectrum indicate that the expan- 
sion of our present universe is accelerating rather than slowing down [1] . This cosmic acceleration 
can not be explained by the four known fundamental interactions in the standard models, which 
is the greatest challenge today in the modern physics. The proposals that have been put forward 
to explain this observed phenomenon can basically be classified into two categories. One is to 
assume that in the framework of Einstein's general relativity (GR), an exotic component with 
negative pressure called dark energy (DE) is necessary to explain this observed phenomena. For 
a good review on the dynamics of different DE models see [2] and references therein. Another 
alternative to account for the current accelerating cosmic expansion is to modify GR theory. 
The well-known modified gravity theories are, for examples, f(R) theory, scalar-tensor theory 
(including Brans-Dicke theory), braneworld scenarios (such as DGP, RSI and RSII), f{Q) theory 
(Q is the Gauss-Bonett term), Hofava-Lifshitz theory, MOdified Newtonian Dynamics (MOND), 
and so forth. For some relevant reviews sec [3]. 

Recently, a new modified gravity theory namely the so-called f(T) theory [4]- [13] was pro- 
posed to describe the present accelerating expansion of the universe without resorting to DE. It 
is a generalization of the teleparallel gravity (TG) [14] by replacing the so-called torsion scalar 
T with f(T). TG was originally developed by Einstein in an attempt of unifying gravity and 
electromagnetism. The basic variables in TG are tetrad fields e^, where the Weitzenbock con- 
nection rather than the Levi-Civita connection was used to define the covariant derivative. As 
a result, the spacetime has no curvature but contains torsion. A vector in TG is parallel 
transported along a curve if its projection Vi = e-i^V^ remains constant, this is the so-called 
teleparallelism. The main advantage of f(T) theory is that its field equations are the second 
order which are remarkably simpler than the fourth order equations of f(R) theory [8]. 

The other interesting issue in modern cosmology is the thermodynamical description of the 
accelerating universe driven by DE or modified gravity. It was shown that by applying the 
first law of thermodynamics (Clausius relation) — dE = Ta&Sa to the apparent horizon ta, the 
Friedmann equation in the Einstein gravity can be derived if we take the Hawking temperature 
Ta = l/(27rfA) and the entropy Sa = on the apparent horizon, where A is the area of the 
horizon [15]. Here, dE is the amount of energy flow through the fixed apparent horizon. The 
equivalence between the first law of thermodynamics and the Friedmann equation was also found 
for gravity with Gauss-Bonnet term, the Lovelock gravity theory and the braneworld scenarios 
[15, 16, 17]. 

Note that in the thermodynamics of the apparent horizon in the Einstein gravity, the geo- 
metric entropy is assumed to be proportional to its horizon area, Sa = [15]. However, this 
definition is changed for other modified gravity theories. For instance, the geometric entropy in 
f(R) gravity is given by Sa = 4cf [-^]> wnere the subscript R denotes a derivative with respect 
to the curvature scalar R. In f(T) gravity, it was shown that when Jtt is small, the first law of 
black hole thermodynamics is satisfied approximatively and the entropy of horizon is Sa = 4cf 
[19], where the subscript T denotes a derivative with respect to the torsion scalar T. 

Besides examining the validity of the thermodynamical interpretation of gravity by expressing 
the gravitational field equations into the first law of thermodynamics in different spacetimes, 
it is also of great interest to investigate the validity of the generalized second law (GSL) of 
thermodynamics in the accelerating universe [20]- [28]. The GSL of thermodynamics like the 
first law is an accepted principle in physics. 

Here, our aim is to investigate the GSL of thermodynamics in the framework of f(T) gravity 
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for a spatially flat Friedmann-Robertson- Walker (FRW) universe filled with the pressureless 
matter. To do this, in section 2, we briefly review the f(T) gravity. In section 3, we investigate 
the GSL of thermodynamics on the dynamical apparent horizon with the Hawking temperature. 
In section 4, for two viable f(T) models we first calculate the effective equation of state and 
deceleration parameters. Then, the null and strong energy conditions are investigated. Also the 
viability of both models is checked by cosmography. Finally, the validity of the GSL is examined. 
Section 5 is devoted to conclusions. 



2 Brief review of f(T) gravity 

In the framework of f(T) theory, the action of modified TG is given by [7] 

'=i6b/ d M'< r > +i 4 w 

where e = det(e* ) = y/—g and L m is the Lagrangian density of the matter inside the universe. 
Also is the vierbein field which is used as a dynamical object in TG. 

The modified Friedmann equations in the framework of f(T) gravity in the spatially flat 
FRW universe are given by [8, 10] 

H 2 = —( Pm + p T ), (2) 
H + \h 2 = -4irG( Pm + PT ), (3) 

where 

^ = ii^ 2T/ r-/-n (4) 

PT = -jtG [ ~^ TfTT + (2T " Ati)h ~ f + ^ ~ T] ' (5) 
T = -6H 2 , (6) 

and H = a /a is the Hubble parameter. Here p m and p m are the energy density and pressure 
of the matter, respectively. Also px and pt are the torsion contributions to the energy density 
and pressure. The energy conservation laws are still given by 

p m + 3H(p m + p m ) = 0, (7) 

Pt + 3H(pt+Pt) = 0- (8) 

Note that if f(T) = T, from Eqs. (4) and (5) we have p T = and p T = then Eqs. (2) and (3) 
transform to the usual Friedmann equations in GR. 

The effective equation of state (EoS) parameter due to the torsion contribution is defined as 
[8, 10] 

Pt 1 T ( 2Tf TT + f T -l \ 

Note that for the de Sitter universe, i.e. H = = T, we have cot = — 1 behaving like the 
cosmological constant. 

Using Eqs. (2), (4) and (6) we have 

(/-2T/ T ). (10) 
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Here, we consider a spatially flat FRW universe filled with the pressureless matter, i.e. p m = 0. 
We ignore the contribution of radiation which is compatible with the observations [29]. From 
Eqs. (2) to (5) for p m = 0, one can obtain 

JT + 2TJtt 

Substituting Eq. (10) into (11) and using T = —12HH, one can get 

T = m( !~ 2T Jj ). (12) 



. /t + 2T/tt 

With the help of above relation, the effective EoS parameter (9) can be rewritten as 

u, f/T-f T + 2Tf TT 

T (/ T + 2T/ TT )(//T-2/ T + l)- 1 > 

For the deceleration parameter 

9 = " 1 -^. ( 14 ) 

using Eqs. (6) and (12) one can obtain 

For the case f(T) = T, we get q = 0.5 corresponding to the matter dominated phase. 

3 GSL in f(T) gravity 

Here, we investigate the validity of the GSL for a spatially flat FRW universe filled with the 
pressureless barionc matter (BM) and dark matter (DM). According to the GSL, the entropy 
of matter inside the horizon plus the entropy of the horizon must not decrease in time [15]. We 
assume that the boundary of the universe to be enclosed by the dynamical apparent horizon f a- 
Also the Hawking temperature on the apparent horizon fx is given by [15] 

T ^ik( l -4k)< (16) 

where 2H? A < ^ ensures that the temperature is positive. 

The entropy of the universe including BM and DM inside the dynamical apparent horizon 
is given by Gibb's equation [20] 

T A dS m = dE m + Pm dV, (17) 

where p m = and V = Airf\ /3 is the volume containing the pressureless matter with the radius 
of the dynamical apparent horizon fx- Also 

E m = — g — Pm- (18) 

Taking time derivative of both sides of Eq. (18) and using (7), Gibb's equation (17) yields 

T A S m = ^lp m (f A -m A ), (19) 
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where the dot denotes a time derivative. Substituting Eq. (10) into (19) gives 



T A S m = ^|(f a - Hf A )(f - 2Tf T ). 



(20) 



Now we need to calculate the contribution of the apparent horizon entropy. Following [19], when 
fxr is small, the horizon entropy in f(T) gravity is given by 



Sa 



Ah 
AG ' 



(21) 



where A = 4irr A . For the special case f(T) = T, Eq. (21) recovers the Bekenstein-Hawking 
entropy Sa = in the Einstein gravity. 

Using Eq. (16), the evolution of horizon entropy (21) is obtained as 



T A S A = t^1 



2G 



2Hr A , 



(2f a/t + r A T/ TT ) . 



(22) 



Inserting Eq. (12) into (22) gives 



T A S A 



1 

2G 



1 



TA 
2Hr A 



(2fA ~ 3m A )2Tf T f TT + 3Hr A f(T)f T T + 2f A /f 
/t + 2T/tt 



(23) 



For the flat FRW metric, the dynamical apparent horizon is same as the Hubble horizon given 
by [30] 

fA = ^- (24) 
Taking time derivative of Eq. (24), using (6), (10) and (11) one can get 



^A = -: 



/ " 2T/ T 



2T \f T + 2Tf TT/ 

Replacing Eqs. (24)- (25) into (20) and (23) and using (6), one can obtain 



TaStu — 



f ~ 2Tf T 
AGT 2 \f T + 2Tf TT 



(4T 2 f TT -4Tf T + 3f), 



TaSa 



f ~ 2Tf T 
4GT 2 V/t + 2T/ tt 



(4T 2 /| t " Tf T + If) (2T 2 /| T -fx- TMtt) 



(f T + 2Tf TT y 



Adding Eqs. (26) and (27) yields the GSL in f(T) gravity as 



T A S 



tot 



_9_ / f-2Tf T 
8G \f T + 2Tf TT 



4f TT + 



f-2Tf T \ ff T + 5Tf TT 



/t + 2T/tt 



T 2 



(25) 

(26) 
(27) 

(28) 



where 5 to t = S m + 5a is the total entropy due to different contributions of the matter and the 
horizon. Note that in the Einstein TG, i.e. f(T) = T, the GSL (28) yields 



?A«Stot = > 0, 



(29) 



which always holds. In what follows, we investigate the validity of the GSL, i.e. T A 5tot > 0; for 
two viable f(T) models introduced in the literature. 
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4 Two viable f(T) models 



Here, we consider two viable f(T) models proposed by [9, 11] to explain the present cosmic 
accelerating expansion of the universe. The first model is a power law 

f(T) = T + W (-T) n , Model 1, (30) 

where \i\ and n are constants [9, 11]. The second model has an exponential dependence on the 
torsion scalar as 

f(T) =T-^ 2 r(l-e /3 ^), Model 2, (31) 

where fi2 and j3 are two model parameters, and To = —GHq [9]. 

The parameters fi\ and /i 2 can be obtained by inserting Eqs. (30) and (31) into the modified 
Friedmann Eq. (2). Solving the resulting equations for the present time gives 

mi = (^t) (e^o 2 ) 1 "", (32) 

1 — Q mo 

& ~ i_(i_2/V ( } 

where Q mn = 87r f^"o i s the dimensionless matter energy density and the index denotes the 
value of a quantity at the present. Note that both models can unify a number of interesting 
extensions of gravity beyond the standard GR. For instance, model 1 for the cases n = and 
n = 1/2 reduces to the ACDM and DGP [31] models, respectively. Model 2 acts like the ACDM 
for f3 = 0. Also both models can satisfy the condition 

lim f/T -> 1, (34) 

T->oo 

at high redshift which is compatible with the primordial nucleosynthesis and CMB constraints 
[9, 11]. For model 1 to be a viable model compared to current data one needs n <C 1 [9, 11]. 
The joint analysis of the astronomical data from SNela+BAO+CMB gives the best fit values 
(n = 0.04lg;l!j, n mo = 0.2721°;°!) for model 1 and (/3 = -0.02l°;^, O mo = 0.272l°;°j*j) for 
model 2 at the 95% confidence level (CL) [9]. 

Using the statefinder geometrical analysis and Om(z) diagnostic method, it was shown that 
both model 1 and model 2 evolve from the standard cold DM (SCDM) to a de Sitter phase [9]. 
Also the effective DE for model 2 with (3 ^ is similar to the ACDM both in the high redshift 
regimes and in the future, while for model 1 with this similarity occurs only in the future 

[9]- 

The evolution of the effective EoS parameter, Eq. (13), for model 1 and model 2 is plotted 
in Figs. 1 and 2, respectively. Figures show that from the early times to the present epoch, both 
models behave like quintessence DE, i.e. lot > — 1, and their effective EoS parameters can not 
cross the phantom divide line [32]. This result has been already obtained by [9, 11]. At early 
times [T/Tq — > +oo), for model 1 and model 2 we have lot — —0.96 and —1, respectively. Model 
2 behaves like the ACDM at high redshift. At the present time {T/Tq = 1), for model 1 and 
model 2 we obtain lot — —0.989 and —0.992, respectively. An interesting result which is absent 
in [9, 11] is that in the future (1 < T/T < 0) for model 1 and model 2 at T/T ~ 0.718 and 
0.719, respectively, we have a transition from the quintessence state, lot > — 1, to the phantom 
regime, lot < — 1- 

The evolution of the deceleration parameter, Eq. (15), for model 1 and model 2 is plotted 
in Figs. 3 and 4, respectively. Figures illustrate that for model 1 and model 2 in the past 
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T/Tq ~ 2.194 and 2.199, respectively, we have a cosmic deceleration q > to acceleration q < 
transition which is compatible with the observations [33]. At early times (T/Tq — > +oo), for 
both models we obtain q ~ 0.5 which indicates that the universe was experiencing a phase of 
deceleration at the early stage of its evolution due to the domination of the matter component. 
At the present time (T/T = 1), for model 1 and model 2 we get qo ~ —0.580 and —0.583, 
respectively, which are in good agreement with the recent observational result — 1.4 < qo < —0.3 
[33]. 

The evolution of px + Pt versus ^- for model 1 and model 2 is plotted in Figs. 9 and 10, 
respectively. Figures show that the null energy condition (NEC), i.e. pt+Pt > 0, is violated for 
model 1 when T/T G (0.025,0.718) and for model 2 when T/T G (0,0.017) and (0.138,0.719). 
When the EoS parameter of DE is less than —1, the universe reaches a Big Rip singularity within 
a finite time. In this case, the NEC is violated. Barrow [35] showed that a different type of future 
singularity, the so-called sudden singularity, can appear at a finite time even when the strong 
energy condition (SEC), i.e. px + 3pr > and pr + Pt > 0, is satisfied. This type of future 
singularity corresponds to the one in which the pressure density pt diverges at T = T s , but the 
energy density pt is finite. The properties of different future singularities in the DE universe 
have been investigated in ample detail by [36]. Our numerical results show that the sudden 
future singularity appears in model 1 and model 2 when (T s /Tq = 0.025, 16irGpT = 20754) 
and (T s /T = 0.138, l6irGp T = 19822), respectively. Also the SEC for model 1 and model 2 is 
satisfied only when T/T G (0,0.025) and T/T G (0.018,0.138), respectively. 



4.1 Cosmographic analysis 

Here, following Capozziello et al. [37] we use a cosmographic analysis to check the viability of 
model 1 and model 2 without the need of explicitly solving the field equations and fitting the 
data. In this approach, the parameters of a given f(T) model must be chosen in such a way that 
the model-independent constraints on the cosmographic parameters (h, qo,jo, sq, lo) obtained by 
fitting to SNela Hubble diagram and BAO data are satisfied. Here, h is the Hubble constant 
and qo, jo, sq, Iq are the deceleration, jerk, snap, and lerk parameters, respectively. Let us first 
start with model 2. Using Eqs. (4.21) and (4.22) in [37] one can find 

1*2= Q-l)(l-O mo ). (35) 

Now the ft = /W(T )/(6 J ff 2 )-( i - 1 ) value for i = (2,3,4,5) where /«(T) = (Pf/dT can be 
expressed as function of (3 only when we fix Cl mo = 0.1329//i 2 from the WMAP7 data. Following 
[37] for each fi value of the sample obtained above from the cosmographic parameters analysis, 
we solve feiP) = H- This yields /3 = —0.2 which takes place in the 95% CL from the model- 
dependent constraints [9]. Then, we estimate the theoretically expected values for the other 
derivatives (/3,/4,/s). The median and 68% and 95% confidence ranges are obtained as 

f — n 071 +0.117+0.148 
/3 — U.O<±_ . 351 _o .914 

f — 1 oocr+0.380+0.473 
J4 — !.OC50_ 1 .302-3 .709 

f — a 447+1-498+1.832 /oe\ 

h - o.447_ 6 . 03 3_ 16L852 . (36) 



Now we compare the above results with the model-independent constraints on the fi values 
given in Table II in [37]. Following [37] we use only the 68% CL which we compare the above 
constraints to. Our comparison shows that the values of (^3, /4, /s) take place in the 68% CL in 
Table II in [37]. Therefore, we conclude that model 2 is favored by the observational data. 
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For model 1, its cosmographic analysis has been already done by Capozziello et al. [37]. 
They showed that for the best fit value of the cosmographic parameters, solving fain) = f 2 
yields n = —0.011. Note that the value of n takes place in the 95% CL from the model- 
dependent constraints [9]. Capozziello et al. [37] found quite small values for (/3,/4,/s) as 
expected for the ACDM model. Here, the disagreement with the constraints in Table II in [37] 
may be due to this fact that depending on the value of n, assuming /W(T ) = for i > 6 m 
Taylor expanding f(T) can fail for model 1 so that the constraints on fi should not be considered 
reliable [37]. 



4.2 Examining the GSL for model 1 and model 2 

Here, we examine the validity of the GSL, Eq. (28), for both models. First, we need to check 
the validity of the horizon entropy relation (21). Because as we already mentioned, Eq. (21) is 
valid only when /tt is small [19]. To check this, we plot Jtt versus T/Tq for model 1 and model 
2 in Figs. 5 and 6, respectively. Figures show that the Jtt is very small throughout history of 
the universe. This confirms the validity of Eq. (21) for both models. 

Now we can calculate the GSL, Eq. (28), for both models. For model 1, the resulting GSL 

is 



T A S- 



tot 



l + /ii(l-2n)(-T) 



n-l 



\n— 1 



(l + W n(l-2n)(-T) r 
1 + (l - n(2 + n))/ii(-T) n - 1 + n(l - 2n)(i - n(9 - 4n))^ 2 (-T) 2n - 2 ] ,(37) 



and for model 2 we obtain 
7A<Stot = - 



8T 



(1 - M2 )T 2 + MT 2 - /3T T + 2/3 2 T 2 )e- 



2 ' 



(38) 



where 



9(1 - mTT* + 9^ 2 [(3 - 2fi 2 )T° - 5/3T (l + fi 2 - fi z 2 )T q 
+/3 2 T 2 (1 + 17 H2 ~ 9/xl)T 3 - 4^ 2 /3 3 T 3 ((5 - 2fi 2 )T 

+^r (i-2/x 2 ))r + i6/x|r 5 ^ 



£2ii 
e t 



(39) 



The evolutions of the GSL (37) and (38) are plotted in Figs. 7 and 8, respectively. Figures clear 
that for both model 1 and model 2, the GSL is satisfied from the early times to the present 
epoch. At early times {T/Tq — > +oo), for both models we have GT\Stot — 1.125. At the present 
time (T/Tq = 1), for model 1 and model 2 we get GTAS to t — 0.116 and 0.010, respectively. 
Note that in the future (1 < T/Tq < 0), Fig. 7 clears that the GSL for model 1 is violated, 
i.e. GT A 5tot < 0, for T/T G (0,0.133) and (0.569,0.718). Also in the future, according to Fig. 
8 the GSL for model 2, like model 1, is violated for T/T G (0.005,0.017), (0.039,0.321) and 
(0.603,0.719). Figures 7 and 8 don't show the negative values of the GSL, because the axis of 
GTA^ot is based on a logarithmic scale. 



5 Conclusions 

Here, we studied the GSL in the framework of f(T) gravity. Among other approaches related 
with a variety of DE models, a very promising approach to DE is related with the modified 
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TG known as f(T) gravity, in which DE emerges from the modification of torsion. The class of 
f(T) gravity theories is an intriguing generalization of Einstein's new GR, taking a curvature- 
free approach and using a connection with torsion. It is analogous to the f(R) extension of 
the Einstein-Hilbert action of standard GR, but has the advantage of the second order field 
equations [11]. We investigated the GSL on the Hubble horizon with the Hawking temperature 
for a spatially flat FRW universe filled with the pressureless matter. For two viable f(T) models 

containing f(T) = T + m(-T) n and f(T) = T — ^ 2 T(1 - e^), we first calculated the effective 
EoS and deceleration parameters. Interestingly enough, we found that for both models there 
is a transition from the quintessence state, ut > —1, to the phantom regime, < — 1, in the 
future. Also both models showed a cosmic deceleration q > to acceleration q < transition 
in the near past which is compatible with the observations [33]. Furthermore, we investigated 
the NEC and SEC and concluded that both models show a sudden future singularity. Also our 
cosmographic analysis cleared that model 1 is unable to predict the observationally motivated 
if 2, /3, /4; h) values but model 2 is favored by the observational data. Finally, we examined the 
validity of the GSL for the selected f(T) models. We concluded that the GSL is satisfied for 
both models from the early times to the present epoch. But in the future, the GSL is violated for 
model 1 when T/T G (0,0.133) and (0.569,0.718) and for model 2 when T/T G (0.005,0.017), 
(0.039,0.321) and (0.603,0.719). 
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Model 1 (n = 0.04, p x = -17240) 
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Figure 1: The evolution of the effective EoS parameter, Eq. (13), versus ^ for model 1, Eq. 
(30). Auxiliary parameters are: n = 0.04, tt mo = 0.272 [9], H = 74.2 Km S" 1 Mpc" 1 [34]. For 
these values one finds m = (^r) {§ H l) l ~ n = -17240. 
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Figure 2: Same as Fig. 1, for model 2, Eq. (31). Auxiliary parameters are: f3 



9. 



0.272 [9], H = 74.2 Km S^ 1 Mpc" 1 [34]. This values gives /j 2 = ^^yn 



mo 



= -0.02, 
-37.5. 
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Model 1 (n = 0.04, m = -17240) 




Figure 3: The evolution of the deceleration parameter, Eq. (15), versus ^ for model 1. Auxiliary 
parameters as in Fig. 1. 
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Figure 4: Same as Fig. 3, for model 2. Auxiliary parameters as in Fig. 2. 
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Model 1 (n = 0.04, m = -17240) 
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Figure 5: /tt versus £- for model 1. Auxiliary parameters as in Fig. 1. 
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Figure 6: Same as Fig. 5, for model 2. Auxiliary parameters as in Fig. 2. 
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Figure 7: The evolution of the GSL, Eq. (37), versus for model 1. Auxiliary parameters as 
in Fig. 1. 
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Figure 8: The evolution of the GSL, Eq. (38), versus tjt for model 2. Auxiliary parameters as 
in Fig. 2. 
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Figure 9: The evolution of pr + Pt versus ^ for model 1. Auxiliary parameters as in Fig. 1. 
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Figure 10: Same as Fig. 9, for model 2. Auxiliary parameters as in Fig. 2. 
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